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ABSTRACT 

Polarization has been detected at early times for all types of supernova, indicating that such systems 
result from or quickly develop some form of asymmetry. In addition, the detection of strong line po- 
larization in supernovae is suggestive of chemical inhomogeneities ("clumps") in the layers above the 
photosphere, which may reflect hydrodynamical instabilities during the explosion. We have developed 
a fast, flexible, approximate semi-analytic code for modeling polarized line radiative transfer within 3-D 
inhomogeneous rapidly-expanding atmospheres. Given a range of model parameters, the code gener- 
ates random sets of clumps in the expanding ejecta and calculates the emergent line profile and Stokes 
parameters for each configuration. The ensemble of these configurations represents both the effects of 
various host geometries and of different viewing angles. We present results for the first part of our survey 
of model geometries, specifically the effects of the number and size of clumps (and the related effect 
of filling factor) on the emergent spectrum and Stokes parameters. Our simulations show that random 
dumpiness can produce fine polarization in the range observed in SNe la (~l-2%), as well as the Q-U 
loops that are frequently seen in all SNe. We have also developed a method to connect the results of our 
simulations to robust observational parameters such as maximum polarization and polarized equivalent 
width in the line. Our models, in connection with spectropolarimetric observations, can constrain the 3-D 
structure of supernova ejecta and offer important insight into the SN explosion physics and the nature of 
their progenitor systems. 

Subject headings: supernovae; spectropolarimetry; numerical modeling 



1. Introduction 

Our understanding of supernovae (SNe) has an im- 
pact on a wide variety of astronomical fields, from star 
formation, interstellar medium dynamics and chemi- 
cal evolution of galaxies through cosmological mod- 
els and the distance ladder. Indeed, our dependence 
on the uniformity of Type la SNe (SNe la) requires 
a greater precision of theory than is generally neces- 
sary in astronomy. Yet there remain a number of unan- 
swered questions in the field that may impact the ac- 
curacy of results obtained from our current models of 
supernovae. One important aspect of SNe not yet well 
understood is their asymmetry, which may impact their 
observed brightness as well as the deposition of kinetic 



energy and enriched material into their environment. 

We have broad indications of SN asymmetry both 
from theory and from observation. High proper- 
motion pulsars from SNe, 7-ray bursts and asymmet- 
ric remnants s upport the idea observationally for core- 
collapse SNe ( Lvne & Lorimenl994 ; Wooslev & BloomI 
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theories of SNe la pro genitors imply breaking o f sym- 
metry due to binarity dWang & Wheelerll2008l) . as do 
the hydrodynamical instabilities and and off-center 
ignition conditions seen in explosion simulations 
( Khokhl ovetalJI 19991; IPlewa et alj|2004l; iRopke et al. 
12007 ). Further, basic expectations regarding rotation 
and magnetic fields suggest asymmetries in the pro- 
genitor star structure with potentially complicated ef- 
fects on the resulting ejecta. 

Observational probes of SN asymmetries have been 
challenging. All SNe yet seen in the modern era of 
astronomical detection have not been resolvable for 
years after the explosion - the first resolved images 



of SN 1987 A, the closest in modern t imes, with the 
Hubble Space Telescope came in 1991 (IPodsiadlowski 
I199I). Further, the most common diagnostic of asym- 
metry in unresolved sources - line profiles - are not 
practical in SNe, where ejecta velocities are typically 
on order of 10,000-20,000 km s^^ or more. Therefore 
the most effective way of probing the asymmetries in 
SNe is with spectropolarimetry. 

Polarization tells us about the asymmetries of an 
unresolved object. In SNe, polarization is generally 
thought to arise from electron scattering in the last few 
optical depths within the photosphere. If the source 
is completely symmetric, all polarization angles will 
be present in equal quantities, and we will measure no 
net polarization - the equivalent of "natural light." If, 
on the other hand, the source photosphere is asym- 
metric, or if some of the polarized light is blocked 
by an asymmetric opacity above the photosphere, dif- 
ferent polarization angles will be present in differ- 
ent streng ths, and we will detect a net polarization 
Kasen et al. (2003) . 



5/7ecfropolarimetry has the potential to give us even 
more insight into the structure of SN. In supernovae, 
the velocity of individual mass elements are set within 
the first few days after the initial explosion. The mass 
of the ejecta quickly stratifies into a pseudo-Hubble 
flow, where distance from the center of the explosion 
is proportional to velocity. Thus, in a SN spectrum, ob- 
servations at a given wavelength will correspond to the 
same mass-elements overtime. By measuring the vari- 
ation of polarization with wavelength, we can observe 
asymmetries as a function of radius, as well as the ef- 
fects of depolarization by line scattering, giving us in- 
formation about chemical and density inhomogeneities 
within the ejecta. 

Many observers have contributed to our knowl- 
edge of t he spectropolarimetric characteristics of 
SNe (e.g., ' """ ' "■' 



Leonard et al. 2006; Maund et al. 
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The results are well discussed in Wang & Wheele 



(l2008h . so here we will only state the most general 
conclusions: Core-collapse supernovae show evidence 
of an underlying axisymmetric explosion mechanism, 
and also show significant deviations from that axisym- 
metry in individual cases. SNe la show little contin- 
uum polarization and thus are likely close to globally 
symmetric. SNe la absorption lines, on the other hand, 
are often associated with strong polarization peaks. 
This can be interpreted as being due to clumpy struc- 
tures in the outer layers, composed of intermediate 



mass elements, which asymmetrically block the pho- 
tosphere. 

To understand the host geometry that produced a 
given spectropolarimetric signature, it is necessary 
to model radiative transfer through possible ejecta 
structures and match the results to observations. The 
problem of modeling SN spectropolarimetry has been 
approached in a variety of ways over the past few 
decades. Most codes for modeling spectropolarimetric 
radiative t ransfer within the ejecta of SN e la, such as 
HYDRA JHoflichlbOOJ I). and SEDONA / Kasen et al. 



2006) use a Monte Carlo approach to line and contin- 
uum radiative transfer and 7-ray deposition within as- 
pherical ejecta to calculate time-dependent spectropo- 
larimetry from a variety of viewing angles. Another 
approach to SN modeling, developed by Jennifer Hoff- 
man, specifically addresses the interaction of core- 
colla pse SNe with circum stellar material in Type lln 
SNe (lHoffmanetaDl2005h . 

These codes are highly detailed and powerful. Yet 
this complexity also means that they are computation- 
ally intensive, making large-scale parameter studies 
impractical. Since polarization is a second order ef- 
fect of the geometry, different source configurations 
may produce the same polarization signature. It is 
therefore necessary to understand the statistical like- 
lihood of different source configurations producing a 
given observational signal. This requires systematic 
modeling of the expected signal given a wide range 
of source geometry. The large number of model con- 
figurations required for reasonable statistics, the three- 
dimensional nature of the problem, and the wavelength 
dependent polarization effects combine to pose a sig- 
nificant computational challenge. 

Since our goal is to perform just such a parameter 
study, we have taken a different approach to supernova 
spectropolarimetric radiative transfer modeling. In the 
interest of simplicity, we have chosen to focus our sim- 
ulations on one potential source of SN asymmetry: re- 
gions of enhanced line opacity ("clumps") in the layers 
above the SN photosphere suggested by the detection 
of strong line polarization in SNe. If present, the na- 
ture of this dumpiness has implications for explosion 
models of SNe. 

The goals of this paper are therefore to calculate the 
line polarization from inhomogeneous distributions of 
elements in the ejecta, and to statistically examine how 
it relates to clump properties. We therefore systemat- 
ically study the effects of various numbers and sizes 
of clumps, placed randomly to recapitulate the effects 



of a range of explosion scenarios as well as those of 
Hne-of-sight orientation. To do this, we have devel- 
oped a fast, flexible approximate semi-analytic code 
for modeling radiative transfer within 3-D inhomoge- 
neous ejecta. In this paper, we present this code and 
our methods for statistically analyzing large numbers 
of simulated spectra, and use them to explore a variety 
of possible clumpy ejecta scenarios. From the ensem- 
ble of resultant Stokes spectra we can determine what 
ranges of parameters predict observable characteristics 
that can constrain aspects of the host geometry. 

In we describe our code and how we extract ob- 
servable trends from an ensemble of calculated spec- 
tra. In §3, we present the results of our analysis of 
the effect of number of clumps and clump size. In §4 
we summarize the current work and discuss the code's 
future potential. 

2. Radiative Transfer Code and Analysis Methods 
2.1. Model Geometry and Assumptions 

The most basic simplifying assumption in the mod- 
eling of supernova ejecta is that all the material is mov- 
ing purely radially and with the velocity of the mate- 
rial increasing linearly with distance from the center of 
the explosion. We can make this assumption because 
the energy that feeds the SN expansion is extremely 
large, and most is input almost instantaneously (on or- 
der of a few seconds). Initial shocks and entrained ra- 
dioactive material will affect the structure of the ejecta 
for a few days, but such effects will be minimal by 
the time the S N is detected, on order of 10 days after 
the explosion (lArnetlill996l) . Further, the expanding 
ejecta is unlikely to encounter any significant circum- 
stellar material at this point in expansion, as that ma- 
terial appears to be cleared out by radiation and other 
pre-explosion dynamics in the years or centuries be- 
fore the explosion. The main exception to this may be 
SN Iln (IHoffman et al.ll200& . and references therein), 
which we do not address in our models. The motion of 
the ejecta can therefore be considered ballistic. 

Our model therefore begins with an N^ grid ex- 
tending from — dmax to +(iijiax in each dimension in 
velocity space (see Fig. [ij. By the time a SN be- 
comes observable, the ejecta will have stratified by ve- 
locity so that this grid is equivalent to a spatial grid at 
a given time. Because of the linear dependence of ve- 
locity on distance, this structure is sometimes referred 
to as a pseudo-Hubble flow, and has many of the same 
characteristics of Hubble expansion. Our supernova is 



centered at d = and has a photosphere at dphot, out- 
side of which is the ejecta through which our radiative 
transfer is calculated. The SN ejecta is first initialized 
to a symmetric, smooth optical depth that decreases 
as a power law from the photosphere (p = po'^"^^ )■ 
(See Table [T] for common parameter values used in all 
simulations reported in this paper.) Regions of en- 
hanced opacity (clumps) are then added, with their 
particular parameters selected randomly from a range 
set for a given simulation. Each simulation consists 
of 1000 spectra in an "ensemble" of realizations with 
the same parameter ranges. Clump parameters that can 
be randomly assigned in this way include the number, 
central optical depth, radii and distance from the pho- 
tosphere of the clumps (rici, Td, Vd and dd, respec- 
tively). We do not specifically vary line-of-sight in our 
results, because given a large ensemble, the random 
placement of clumps will recapitulate a variety of lines 
of sight. This simplification greatly reduces computa- 
tional and analytical complexity. 



Nsrid 




Fig. 1 . — Model grid geometry. 



Our Stokes spectra have N wavelength bins cen- 
tered around a central wavelength, in this case 8000 A. 
Note that the exact wavelength is only relevant to the 
background blackbody spectrum; for these simplified 
models the "line" is pure trapped resonance scattering 
and particular line physics are not addressed. N here 
is the same as the grid dimension. 

Our radiative transfer is done in two parts. We begin 
with the intensity and polarization of the continuum 
flux, which for a spherically symmetric photosphere 
can be represented in polar coordinates as 



lir,^) \ ( Ip{r) 

Q{r,<P) - /p (r)Pp(r) cos 20 
U{r,4)) ) \ Ip{r)Pp{r)sm2(f> 



(1) 



where /p(r) is the specific intensity and Pp{r) the per- 
cent polarization emerging from an impact parameter 
r. The Stokes parameter corresponding to circular po- 
larization, V{r, (j)), is taken to be zero, and thus not 
included in our equations. This is consistent with ob- 
servations, as no circul ar polarization has yet be en de- 
tected in a supernova (IWang & Wheeler! l2008h . We 
also use use the Stokes parameter notation given in that 
paper, where /, Q and If represent the measured vec- 
tor components, while Q = Q/I and U ~ U / 1 have 
been normalized. Integrating the intensity over a plane 
perpendicular to the line of sight gives the Stokes flux 
at a specific wavelength 
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(2) 

where r^ is the photospheric radius, dmin is the impact 
parameter at which the given plane intersects the pho- 
tosphere (or zero if the photosphere is not intersected) 
and S is the line source function, assumed to be pure 
scattering, S = W{r, (p)Ip where W is the geometrical 
dilution factor. The plane located at distance z along 
the line of sight (with a projected velocity v^) corre- 
sponds to a wavelength A = Ao(l + Vz/c). In theory, 
each of these function could have wavelength depen- 
dence, but for simplicity over such a small wavelength 
range we do not vary any of these functions except by 
letting /p(A) vary according to a blackbody spectrum. 

In our simulations, Ip{r) and Pp{r) are set by the 
results of a spherical Monte Carlo calculation. Note 
that for a spherical photosphere, the polarization at 
the limb is higher than i n the plane parallel calcu- 
lation of [C handrasekhar (196(J), used before in e.g., 
Shapiro & Sutherland ( 1982). The values of the Q and 
U Stokes vectors across the observer-oriented face of 
the photos phere are taken from the model calculated in 
Kasen et a l. (2003). 

Equation |2] represents the effect of line opacity 
above the photosphere in obscuring or diluting contin- 
uum flux is calculated. Operationally, for the radiative 
transfer through the ejecta, we calculate the interacting 
wavelength for each cell in the 3-D velocity grid us- 
ing the Sobolev approximation (Sobolev 1960) for line 
opacity. The light from the photosphere that passes 



through that cell is then attenuated or scattered at the 
resonant frequency by the amount determined by the 
cell optical depth. The combined effects of each cell 
on all wavelengths gives us our line profile, and asym- 
metric depolarization of the photospheric flux will give 
us a spectropolarimetric signal. 

Note that we make two main assumptions in this 
model: First, that there is a clear division between the 
electron scattering photosphere and the line forming 
region of clumps. In reality, there is no sharp divi- 
sion, and some electron scattering will be coincident 
with the line scattering. However, if the clumps can 
be represented as sufficiently detached from the pho- 
tosphere, this model will still provide relevant results 
')S' (see discussion in Kasen et al. 2003) . 

Second, we make the common assumption that all 
line scattered light is depolarized. Though this is cer- 
tainly a simplification of the physical picture, it is a 
reasonable a pproximation in most cases. The reason- 
ing given in iHoeflich et al.l (1199 6) is that collisional 
timescales are less than those for absorption and re- 
emission in the ejecta; thus the random scattering pro- 
cesses will erase incoming geometrical information for 
an absorbed photon. We add that even if collisions are 
not dominant, line polarization is unlikely to be sig- 
nificant becaus e lines are at b est weakly polarizing, 
comparatively (lHamiltonlll947h with the greatest po- 
larization in the side scattering case. In our geometry, 
photons that are side scattered into the line of sight will 
be seen in the emission portion of the line, while ob- 
served polarization is in the absorbed part of the line. 
Further, any light line-scattered into the line of sight 
will be a small contribution to the total flux, making 
such a signal harder to detect. 

Thus in our model, polarized rays of light from 
a pre-computed spherical photosphere are obscured 
by clumps of depolarizing line opacity. Because the 
clumps block the photosphere in an asymmetrical way 
at different velocities, they give rise to polarization 
over the line flux absorption feature. 

Results from our code have been comp ared to those 



of Kasen's previous semi-analytic codes ( Kasen et al. 



2003) with the same host configurations, to good 
agreement. See Hold (120091) for more details. Flux 
and polarization spectra and Q-U diagrams for three 
sample realizations are shown in Figures |3] Sand |5] 



2.1.1. Regimes in nd-rd Parameter Space 

The code can be run with any number or size 
of clumps, but it is worth examining what is physi- 
cally implied by different regions of rtci-7'ci parame- 
ter space. In the small-clump limit, we are modeling 
small regions of enhanced concentration of an ele- 
ment, surrounded by smooth ejecta which contribute 
the bulk of the line profile but no net polarization. 
As the clumps become larger or more numerous, the 
model more closely resembles large scale variations in 
the chemical composition of the overall ejecta struc- 
ture. Eventually, however, the simulations with the 
most and largest clumps will completely fill the model 
grid, and represent a substantial change to the nature of 
the ejecta, rather than a perturbation of it, to a degree 
not supported by observation. 



Thomas et al.l (120021) attempted to constrain clump 



size and filling factor using line flux profiles only. 
Comparing results of a numerical clump model to 
observations of the Si II line, they concluded that 
clumps larger than 0.08 times the photospheric radius 
(1600 km s^^ in our model) would create line-of-sight 
variations larger than what was seen. Their model was 
substantively different from ours, however, as they rep- 
resented clumps as optically thin and surrounded by 
zero line opacity. The Une profile is thus due entirely 
to the clumps, rather than clumps being a perturbation 
on the line structure. Therefore in our model, the line- 
of-sight effects will be smaller even for large clumps, 
and should have less variation with angle as long as 
there are enough clumps to fill a majority of the ejecta. 
In an attempt to quantify the useful parameter limits 
in our model, we use a metric based o n the ph otodisk 
covering factor (PCF) used by .Thomas et al.l (12002) . 



Table 1 : Model parameters used for simulations in this 
paper. Sizes are given in velocities, which are propor- 
tional to distances in the homologous flow of SN ejecta 
(see[T]). 



Parameter 


Value 


N 


100 


^inax 


20,000 km s-i 


f^phot 


10,000 km s-i 


Tbb 


13,000 K 


Tel 


10 


^cent 


8000A 


nej 


8 



which they define as the ratio of the area obscured by 
clumps contained in the projected photodisk area to the 
total photodisk area for each plane of common velocity 
along the line-of-sight. The sum of the ratios over all 
velocities gives a measure of volume filling factor for 
the ejecta. For these simulations, there are 50 velocity 
planes on the observer side of the grid, so the sum of 
the ratios has a maximum of 50 if each plane has 1 00% 
of the available photodisk covered. 

We therefore delineate four regions of parameter 
space in our model, as shown in Fig. |2l 

I The sm all-clump reg i me, r d < 1600 km s~^, 
used bv .Thomas et alj (120020 



II Medium-clump regime with smaller filling factor 
(PFC < 38) 

III Medium-clump regime with larger filUng factor 
(38 <PCF < 48) 

IV Full ejecta (PFC > 48) 

Region I has small clumps and is numerically and 
physically plausible. Regions II and III have more sub- 
stantial dumpiness. In Region II the ejecta is less full, 
which might lead to greater line-of-sight variation than 
is seen in observations. Region III has a higher fill- 
ing factor and thus fewer line-of-sight effects, but also 
results in more clump overlap, implying the possibil- 
ity of more complex polarization effects that may not 
be completely captured with our simplified radiative 
transfer. Region IV represents a full numerical grid 
and is not likely to occur in SNe but is shown here 
to demonstrate the behavior of the model in extreme 
cases. 

2.2. Analysis of Spectra 

Given that each simulation creates an ensemble of 
1000 individual Stokes spectra, we developed tools to 
analyze and visualize the results in a statistical way. 
We therefore calculate certain salient, observable num- 
bers from each spectrum, and determine the "expected 
value" for each configuration. The numbers calcu- 
lated from each spectrum are minimum flux and the 
velocity and polarization where it occurs; maximum 
percent polarization (P, see eqn. O at any frequency 
in the spectrum (Pmax); polarized equivalent width, 
Q- vector equivalent width, U- vector equivalent width, 
and Stokes equivalent width (PEW, QEW, UEW and 
SEW, see eqn. lU; and theoretical covering fraction 
(TCF, eqn.Ell. 
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Fig. 2. — Regions of rd-nd parameter space. Re- 
gion I represents the small-clump limit, with rd < 
1600 km s^^. Regions II and III represent medium- 
sized clumps, with a smaller ejecta filling factor in Re- 
gion II and greater clump overlap in Region III. Region 
IV is the filled-ejecta, theoretical limit of the model. 



The polarization is defined as 



P 



C/2 



Vo^+lP 



(3) 



We also define a polarized version of line equivalent 
width: 



[X]EW = Yl ^"^ 



. , '-^4^^ 1 AA, 



I, 



c.X 



(4) 



where [X] can be P, Q, or U and Ic is the continuum 
flux. Since QEW and UEW should each be sampling 
the same underlying distribution even if their values 
differ in a given spectrum, we add these separate distri- 
butions together to double the sample size. We call this 
combined value "Stokes equivalent width," or SEW. 
In this paper, we focus primarily on the analysis of 
Pniax and PEW results. 

The concept of the "theoretical" covering fraction 
is also worthy of more discussion. It is useful to have 
a way of calculating a covering fraction in order to fa- 
cilitate comparison to the results of explosion codes; 
however the correct measure of covering fraction in 
our simulations is not as simple as one might expect. 
First, our randomly placed clumps may not be in the 
Une of sight, and therefore may not in that sense be 
"covering" the photosphere from the observer's point 
of view. However, to most easily compare our results 
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Fig. 3. — A calculated line profile, spectropolar- 
ization profile and Q-U diagram for the nci=32, 
rci=1600 km s^^ case. This particular realization has 
several clumps in the line-of-sight and therefore has a 
number of broad polarization peaks and distinct varia- 
tions from a single dominant axis in the Q-U plane. 
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Fig. 4. — A calculated line profile, spectropolarization 
profile and Q-U diagram for the same clump configu- 
ration as in figure |3] but with rci=400 km s^^ instead 
of 1600 km s^^. Note that the polarization peaks are 
significantly narrower and > 10 times smaller in mag- 
nitude. 
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Fig. 5. — Another calculated line profile, spectropo- 
larization profile and Q-U diagram for the nci=32, 
rci=1600 km s^^ case (e.g., it is from the same en- 
semble as the example in figure |3]l. This particu- 
lar realization produces a classic "loop" in the Q-U 
plane. The difference that leads to this loop is hinted 
at by the plots of percent polarization: in figure [3] the 
clumps are separated in velocity space, so the Q-U di- 
agram shows distinct deviations from zero, with a sep- 
arate position angle for each clump. In contrast, the 
clumps in this realization are overlapping in velocity 
and wavelength, so the Q-U diagram shows a smooth 
varying from one clump's position angle to another, 
creating a "loop" structure. 



with the predictions of explosion codes, we need to ac- 
count for the total number of clumps. Further, because 
our clumps are distributed in velocity space, the cov- 
ering fraction will be different at different frequencies 
depending on where the Sobolev line interaction oc- 
curs. The covering fraction at each wavelength in the 
line of sight can be calculated (e.g., the "photodisk cov- 
ering fraction" of section |2.1. Il l however this measure 
is not ideal here because it is a) a spectrum rather than 
a single number, b) subject to line-of-sight effects and 
c) more complicated to compare to explosion models. 
We therefore have chosen to use the theoretical cover- 
ing fraction, which we define as follows: 



TCF = 



t: 4-^ci 



(5) 



In words, we find the fraction each clump subtends 
of the total area of the sphere located at the distance of 
the clump, and then take the sum of these fractions for 
all clumps. There are limitations to this method, par- 
ticularly in that our clumps may overlap either within 
the opacity grid or in line of sight, and may therefore 
overestimate the covering fraction, potentially creat- 
ing a TCP that is greater than one, even if there are 
some locations or wavelengths where the photosphere 
is unobscured. It does, however, create a single, con- 
ceptually clear number to use in calculations and com- 
parisons. 

2.3. Analysis of Ensembles 

Once the appropriate metrics are calculated for each 
individual spectrum in an ensemble, we can perform 
statistical analysis on the results. If we have 1,000 re- 
alizations of a given set of parameters, we have 1,000 
predictions of, for instance, the maximum polariza- 
tion that could be produced by that configuration. A 
Gaussian fit to a histogram of the maximum polar- 
izations thus produced allows us to predict the like- 
lihood of a given polarization being produced by that 
geometry. We can also use such model parameters as 
number of clumps (rid), clump size (rd) or TCP to 
connect these geometries to those predicted by differ- 
ent SNe explosion modeling codes, and thereby deter- 
mine what spectropolarimetric signatures these explo- 
sion codes would imply in our simulations. 

In Pigure|6]we show some of the statistical results 
for an ensemble, namely the histograms of Pmax and 
SPW. In this ensemble, clumps were randomly dis- 
tributed between the photosphere and the edge of the 



grid, and with nci=16 and rci=1600 km s ^. Addi- 
tional parameters are as in Table [T] 

Analysis of such histograms can lead to insights 
into the factors contributing to the overall distribu- 
tion. In this ensemble, we see three main contribu- 
tions. Pirst, the peak at zero polarization represents a 
combination of symmetrically placed clumps and in- 
stances where the clumps were randomly placed so 
that none of them lie in the line-of-sight. Between zero 
and ^0.15% we see a steadily increasing slope that 
make up the bulk of the asymmetric population, fol- 
lowing the expected properties of clump distribution. 
Third, in addition to this bulk population, there are 
a small number of outliers in the distribution. These 
outliers have the potential to demonstrate interesting 
though unlikely configurations, and possibly to distin- 
guish between similar clump models. 

In addition to analysis of individual ensemble dis- 
tributions we seek to understand the way our results 
change with variations of simulation parameters. To 
see trends in metrics between ensembles, we must find 
ways to characterize each distribution in ways easily 
comparable to others with different parameters. One 
measure we use is a Gaussian fit the distribution, which 
gives a center and a width. The true mode of the distri- 
bution is also calculated, and finally, for some metrics, 
we find the mode of the distribution if zero is excluded. 
This second mode can be useful in distributions that 
can have a peak at zero due to, for example, no clump 
falling within the line-of-sight, which is distinct from 
the behavior of the metric when any clumps are in po- 
sition to affect the observed signal. This second, non- 
zero mode would therefore be useful for metrics that 
are not expected to be zero when there are any asym- 
metries present (e.g., polarized equivalent width) but 
not for metrics where expected values should be dis- 
tributed about zero (e.g., Qmax)- 

As an example, we again take the case of Pmax in 
the line as shown in Pigure|6] Here, the Gaussian fit 
is not a complete representation of the underlying dis- 
tribution. Indeed, since P is an intrinsically positive 
quantity, the distribution of P measurements in an en- 
semble should not be Gaussian by definition. Por this 
metric, the mode seems to be a better measure of the 
behavior of the distribution produced by our simula- 
tion. Unfortunately, the mode cannot capture the width 
of the distribution, only its most likely value. 

In contrast, there are some metrics where the distri- 
bution of measurements is more usefully represented 
by the Gaussian fit than by the mode. Both Q and U 



individually are expected to be randomly distributed 
about zero, and thus the combination of the two distri- 
butions should be as well. We therefore expect that this 
data will be more accurately represented by a Gaussian 
fit. In Figure|6]we show that the distribution of Stokes 
equivalent width for the same ensemble is indeed more 
accurately fit by the Gaussian than was the Pmax- Here 
the width of the fit is a more useful quantity than the 
mode, which will be at or very close to zero. 

2.4. Inter-Ensemble Parameter Studies 

Finally, in order to investigate the effect of different 
variables on the predicted spectropolarimetric signals 
we calculate a series of ensembles holding all param- 
eters constant except one. We refer to the resulting set 
of ensembles as a "run." To visualize the effects of that 
parameter, in Figures [8] to [T3] we plot the the Gaus- 
sian fit parameters and modes determined for each 
metric in each ensemble versus the varyied parame- 
ter. The width of the Gaussian is shown by the gray 
area. We plot the true mode for zero-centered metrics 
(i.e. Stokes equivalent width) and non-zero mode for 
metrics that are not expected to be inherently symmet- 
ric. 

3. Results 

We have performed a preliminary survey of the 
effect of three interrelated factors on the spectropo- 
larimetric signatures of SNe: clump size, number of 
clumps and TCF. We note that in addition to likely pa- 
rameter values, for some of these simulations we use 
numbers or sizes of clumps that are at physical ex- 
tremes in order to explore the theoretical limits of the 
models. (For instance, an ensemble with 1024 clumps 
of radius 8000 km s^^ would not be best characterized 
as representing a "clumpy ejecta," but can provide in- 
sight into trends in the model. See Section l2.1.1l ) 

In this paper we present complete results for seven 
runs: two have constant rd and vary rtci (parameters in 
Table |2]l, two have constant rtciand vary rci(Table[3]|, 
and three vary both rdand rtcisystematically to main- 
tain an approximately constant TCF (Table lU. For 
each of these runs we show trends for several metrics. 
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Fig. 6. — The histogram and Gaussian fit for the maxi- 
mum polarization produced at any point in the line pro- 
file (above) and Stokes equivalent width (below) for an 
ensemble with nci=16 and rci=1600 km s~^. In the 
^max plot, the peak at zero represents the combination 
of symmetrically placed clumps and instances where 
no clumps appear in the line of sight. The Gaussian fit 
was set to avoid this peak; the mode of the distribution 
is calculated for the full distribution, and again exclud- 
ing the peak at zero. In this case, both modes are the 
same. Note that the Gaussian fit is a better representa- 
tion of Stokes equivalent width data than for Pmax, as 
expected, because Stokes equivalent width should be 
scattered about zero. 



3.1. Constant r. 



We conducted two runs with constant clump ra- 
dius, one with rci=3200 km s^^ and the other with 
rci=6400 km s^^ (see Table |2] as well as Table [T] for 
complete parameters). 



In our clump model, polarization of SNe is due to 
asymmetric obscuring of the photosphere (see §??). 
We expect an increase in polarization as the number of 
clumps increases. But beyond a certain point, added 
clumps will create symmetries with other clumps, and 
the polarization signal will decrease. Therefore we be- 
gin by looking for this pattern of increasing polariza- 
tion at small nd followed decreasing polarization at 
large rici- While this behavior is not readily apparent 
in the rci=3200 km s^^ run (Figure|7]i, it is seen in the 
rci=6400 km s~^ series (Figure[8]i. This is likely due 
to the fact that larger clumps will cover a broader range 
of Q and U vectors, meaning that it is easier to regain 
symmetry with fewer clumps. 

In the Pmax versus rid plot of Figure |7] we have 
overlaid a fit to the linear portion of the Pmax data, 
which has a slope of 0.33. A similar set of plots for 
rci=6400 km s^ ^ are shown in figure|8l the slope of the 
central approximately linear portion of the Pmax ver- 
sus rici plot is similar, 0.32. 
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3.2. Constant n^ 

We conducted two runs with a constant number of 
clumps, one with rid — 16 and the other with nc\=32 
(see Table [3] as well as Table [T] for complete parame- 
ters). 

Fi gure s l9l and [TOl show that in the constant nd runs 
PmaxOCJ'ci- At larger r^ this relationship breaks down 
due to increasing geometrical effects from clump over- 
lap, both spatially in the line-of-sight and spectrally in 
velocity space. In the Pmax versus rc\ plots of Figures 
|9]and[T0] we have overlaid a fit to the linear portion of 
the data. The slopes are 1.74 and 1.72, respectively. 
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3.3. Constant TCF 

For this series of runs, we held the theoretical cov- 
ering fraction (TCF, see ^ 12. 31 for definition) constant, 
varying r^cifrom 1 to 1024 and r^ proportionately to 
maintain a TCF of « 0.15, « 0.60 and « 2.4. The 
actual TCF obtained for each ensemble is closest to 
the nominal value at higher rid where the value is less 
likely to be skewed by random placements of clumps. 
Figures [TT] [12] and [13] all show that, as anticipated, 
large numbers of small clumps leads to recapturing 
symmetry and a decrease in polarization. 

In Fi gures [TTl [T2l and [T3] we show the log-log plots 
for the TCF- 0.15, TCF ~ 0.60 and TCF ~ 2.4, 
respectively. In these plots, we have overlaid fits to the 
Unear portions of both Pmax plots. In the TCF ~ 0.15 



Fig. 7. — The dependence of Pmax and PEW on rid in 
the rci=3200 km s^^ run. A fit to the linear power- 
law portion (slope is 0.33) is plotted in the Pmax figure 
(above). TCF of the ensemble is shown, at the top of 
each figure. 



10 




TCF 











1 






10 








100 




100.000 
10.000 


; 


0. . 


-.-O ■ 


o- 


..■■o" 






■o. 


o 


0, 




. o ■■■ 


■o . 










s 

LlJ 


1.000 
0.100 
0.010 
0.001 


r 
"1 








f 










- 
1" 



10 100 



1000 



Fig. 8. — The dependence of Pmax and PEW on rid in 
the rci=6400 km s^^ run. A fit to the Hnear power- 
law portion (slope is 0.32) is plotted in the Pmax fig- 
ure (above). Compared to the rci=3200 km s^^ run 
(Figure|7]i we see significant deviation from power-law 
behavior at large numbers of clumps, indicating that 
symmetric areas of the photosphere are being covered 
by the larger clumps. TCF of the ensemble is shown, 
at the top of each figure. 
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Fig. 9. — The dependence of Pmax and PEW on rd in 
the nci=16 run. A fit to the linear power-law portion 
(slope is 1.74) is plotted in the Pmax figure (above). 
TCF of the ensemble is shown, at the top of each fig- 



ure. 
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Fig. 10. — The dependence of Pmax and PEW on r^\ in 
the nci=32 run. A fit to the Unear power-law portion 
(slope is 1.72) is plotted in the Pmax figure (above). 
TCF of the ensemble is shown, at the top of each fig- 
ure. 



Fig. 11. — The dependence of Pmax and PEW on 
ricx and r^\ in the TCF~ 0.15 run. Linear fits are 
oc — 0.66nci and oc 1.32rci. 
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Fig. 12. — The dependence of -Pmax and PEW on 
rici and rd in the TCF~ 0.60 run. Linear fits are 

oc — 0.55nci and oc l.llrd. 



Fig. 13. — The dependence of Pmax and PEW on 
rid and rd in the TCF~ 2.4 run. Linear fits are 

oc — 0.33nci and ex 0.65rci. 
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case, we find that this slope is 1.32 in the rd plot and 
-0.66 in the rid version; for TCF ~ 0.60 they are 1.11 
and -0.55 respectively. In the ~ 2.4 run, a fit to the data 
gives slopes of 0.65 for r^i and -0.33 for nd, though 
the behavior no longer appear to be truly linear, due to 
complicated effects of overlapping of clumps. 

In the constant rd and rid runs ( ij3.1| and ^^, the 
slope of the fit to the Pmax remained the same between 
runs, at least in the portions that are roughly linear 
The TCF runs show less uniformity (even excluding 
the highest TCF run where behavior does not seem lin- 
ear) than in the other cases. This may indicate a trend 
where increased covering fraction leads to a flattening 
of polarization dependence in both rd and rici- 

The difference between the TCF ~ 0.15 and TCF ~ 
0.60 is likely due to small statistical variations at very 
small covering fraction. We therefore take the results 
from the TCF ~ 0.60 run as most representative of the 
regime of power-law behavior for the TCF runs, where 
there are enough clumps for reasonable statistics, but 
not so many that complex overlapping effects make the 
polarization dependence non-linear. 

3.4. Polarization dependence on nd and rd 

In order to understand these results, we seek a sim- 
ple model for the relationship between polarization 
and rid, ?'ci and TCF. We therefore propose the fol- 
lowing simple analytical model and compare it to our 
model's numerical behavior in order to test consis- 
tency, and to provide insight into both simple and com- 
plex regimes. 

We begin our proposed model with the polariza- 
tion due to one clump. Pel- The polarization from that 
clump will depend on the photospheric area covered 
by the clump, ex r^^. But the polarization will also de- 
pend on what portion of the photosphere is covered - 
if the clump is very large, it will eventually cover or- 
thogonal vectors in Q-U space, and therefore reduce 
the polarization due to that clump. Thus we make the 
assumption that the polarization due to the clump will 
be proportional to the area of the clump times an un- 
known function of rd, which we in turn assume to be 
a power of rd: 



Pel ex rlJir^) cc rl,r^, 



(6) 



The polarization should also be a function of rid, 
which should be roughly linear for small covering frac- 
tions. As we increase the number of clumps, we will 
begin to have more than one clump at the same veloc- 



ity, each covering different Q and U vectors. While 
(Q) = 0, (Q^) is not. Like a random walk, Q^ ex rid. 



The same is true for U. Thus, since P = ^/Q^+lP 
(equation|3]l, the total polarization will also increase as 
y^rid. (A similar derivation of the P ex -^nd behav- 
ior is shown in Richardson et al. (1996).) Thus in our 
model, the total polarization will depend on the num- 
ber of clumps and the polarization per clump: 



Ptot oc g{nci)Pc\y^i 



(7) 



Here g{nc\) is an as-yet unspecified function of 
rid representing the non-linear effects on the total 
polarization due to multiple clumps. Assuming that 
(?(rici) oc ri[.[ and combining equation [T] with equation 
|6]this becomes 



Ptot « ni,r^,nlCrl 



i+l/2 ™+2 
"cl 'cl 



(8) 



Equation[8]predicts that within a run, the total polar- 
ization should have a power-law dependence on rd and 
rid. This prediction agrees with the behavior of our 
models shown in i)3.1| and i]3.2| The actual dependence 
can be disentangled by holding either rd or rid con- 
stant: 



Ptot oc 



'cl 

1+1/2 



for const, rid 
for const, r^ 



(9) 



We can now use our simulations to determine the 
values of I and m. The metric that will most directly 
correspond to the Ptot derived above will be Pmax be- 
cause it is less influenced by complicated geometry ef- 
fects in the line. (PEW should also follow this relation- 
ship in cases where polarization is due to a few discrete 
clumps.) 

From our simulations, the slope of P,„ax as a func- 
tion of rici for constant rd is 0.33 (Figure |7]i. Thus 



1+1/2 



,0.33 



, and / 



-0.17. Our results for con- 



stant rid (Figure |9]l show Pmax as a function of rd is 
ex rlC^ = r;'/+^ and so ni = -0.26. 

To check for consistency, we can use our model and 
these results to predict the behavior of Pmax for runs 
with constant TCF. For constant TCF, rid'^d i^ '-^^i- 

Q 1/2 

stant, and ri-d ex r^j and rd ex n^j . Substituting 
these relations into eqn. |8]we find 



Ptot oc 



,1.08 



-0.54 



const. TCF 



(10) 



These predicted slopes agree remarkably well with 



the fit to our numerical results: r^i^^ and ; 
urefTIb. 



-0.55 



(Fig- 
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Note that this impHes that the total polarization is 
decreased by nci~°'^^ if there is substantial clumping. 
Similarly, total polarization is decreased by rd^'^'^^. 
This implies that the larger the radius of the clump, 
the more likely it is to extend into less polarized parts 
of the photosphere or cover complementary Q and U 
vectors, and thereby reduce the resulting polarization 
per clump. 

Similar behavior for the polarization can be derived 
by introducing clumps to the integrals in equations [T] 
and |2] For non-overlapping, optically thick clumps, 

— 1/2 

the average polarization would be propto n^j ' and 
propto rl^. The difference from these values in our 
results is due to the more complicated effects, such as 
overlap and velocity-width of the clumps, that can be 
handled more effectively in the numerical model. 

3.5. Connection to Observations 

One of the greatest obstacles to the effective use of 
spectropolarimetry in astrophysics is the degeneracy 
between source and signal - because polarization is a 
second order effect, there will likely be more than one 
host configuration that will produce a given observa- 
tion. One of the most important goals of this paper is 
to find ways to break polarimetric degeneracy where 
possible, and otherwise to estimate the range of possi- 
ble configurations leading to an observation, and their 
relative likelihoods. 

Two factors that can break degeneracy in polarime- 
try are wavelength and time dependence. For instance, 
while continuum polarization for an ellipsoidal photo- 
sphere may be the same as that for a spheroidal pho- 
tosphere with chemically inhomogeneous ejecta, the 
wavelength dependence of the polarization will be dif- 
ferent. The former case will also likely maintain the 
same degree of polarization and position angle over 
time, while in the latter, individual clumps will contin- 
ually become optically thin over time as the ejecta ex- 
pands, likely resulting in variation in both degree and 
orientation of polarization. 

Our code incorporates the first of these factors in- 
trinsically, namely wavelength dependence in the line 
profile. Degeneracies for a single SN can be further 
constrained using observations at different epochs. 
Take, for instance, the case where one epoch of ob- 
servations indicate that the asymmetry may be due to a 
single large clump covering most of the photosphere, 
or to several overlapping medium-sized clumps at 
about the same distance from the photosphere. As 



the ejecta expands, the outer portions of the clumps 
will become optically thin while the portion of each 
clump closest to the photosphere may remain optically 
thick. In later observations of the single-clump case, 
we would observe the signature of one clump reducing 
in size. With several overlapping clumps, as we see 
further towards the center of the SN the single clump 
would appear to break up and the spectropolarimetric 
signature would become that of several small clumps. 
Our simulations predict that this difference should be 
detected in the relative change in Pmax and PEW at 
different epochs. 



E 




Fig. 14. — The expectation value of Pmax as a func- 
tion of rici and rd predicted by our simulations. Here, 
a measurement of Pmax in a SN line would be a plane, 
constraining the host system's parameters to those de- 
scribed by an intersecting contour in the above plot. 



More directly, our simulations also give us a pow- 
erful method of constraining clump parameters by us- 
ing observational metrics from the same epoch. Using 
our results from across all of our runs, we can create a 
map of a metric, such as Pmax, as a function of rid and 
Td (Figure O. This map will have the form of a to- 
pographical map, though rather than having a single 
value at each point, the surface would have a thickness 
representing the likely range of Pmax from our Gaus- 
sian fit for that (nd,rd). An observational measure- 
ment of Pmax in the line would be a plane, also with a 
thickness representing the uncertainty in the measure- 
ment. That plane would intersect our map, constrain- 
ing the source rici and rd to the contour thus created. 
The intersection of uncertainties from our model and 
observations would give an uncertainty on the predic- 
tion. 
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Fig. 15. — The expectation value of PEW in the line as 
a function of rid and rd predicted by our simulations. 
Here, a measurement of PEW in a SN line would be 
a plane, constraining the host system's parameters to 
those described by an intersecting contour in the above 
plot. See the caption of Figure [14] for a description of 
the dashed and cross-hatched regions. 



This method becomes even more powerful when we 
use more than one metric: the intersection of observed 
polarized equivalent width with our model (FigurefTSt 
would create another contour in r^ and rich and the 
two contours could be combined to further constrain 
source geometry. Each metric would provide another 
constraint and narrow the possible source parameters. 
Consistency between these predictions for a given ob- 
servation also provides a test of how well our models 
represent the physical system of the SN. 

4. Summary 

In this paper, we have presented a new spectropo- 
larimetric radiative transfer code, designed to model 
the effects of chemical inhomogeneities ("clumps") 
in the ejecta of SNe on a pure-scattering line pro- 
file. While most approaches to the modeling of spec- 
tropolarimetry have focused on detailed m odeling of a 
small number of h ost geometries (e.g.. lHoflichll2009t 
iKasen et al.ll2006l) we have instead chosen to use sim- 
plified radiative transfer to model large numbers of 
host configurations (realizations) with a broad range of 
parameters. This approach has two goals: first, to cre- 
ate a catalog of spectropolarimetric line profiles and 
host geometries that can be compared to observations; 
and second to give an estimate of how unique (or de- 



generate) is the connection between model and obser- 
vation. 

This paper presents the results of our first survey 
of the likely polarization signals from clumpy SN 
ejecta, specifically the effects of the number and size of 
clumps (and the related effect of covering fraction) on 
the emergent spectrum and Stokes vectors. We have 
paid special attention to the connection between un- 
derlying clump geometry and predictions for robust 
observational parameters, such as Pmax and polarized 
equivalent width of the line. By themselves, these re- 
sults allow observers to understand the range of possi- 
ble configurations that may produce a given measure- 
ment, including the relative likelihood of each. 

One important result of our simulations is that ran- 
dom dumpiness can regularly produce peak line po- 
larization of 1% with small to medium clumps, and 
up to 1.6% in simulations where we allow clumps to 
largely obscure the photosphere at some wavelengths 
(see 12.1.11 ): polarization produced by individual real- 
izations can produce even higher values. These aver- 
age values are in the range observed for SNe Type la 
(1-2%), and indicate that line polarization is a plausi- 
ble explanation for existing data. 

We have also shown that loops in the Q-U plane, 
which are often obs erved in SNe. can be produced 
by dumpiness alone. iKasen et alJ (120031) showed that 
such loops could be produced by a clump above an as- 
pherical photosphere. Here we have shown that these 
loops can be produced with a spherical photosphere 
but multiple clumps, if clumps with different position 
angles overlap in velocity (see figures [3] and |5]). 

In order to connect our simulations to observations 
more concretely we present the results of our models 
for observable metrics P,nax in the line and PEW as a 
function of 7ici and rd. A measurement of Pmax is a 
flat plane with finite thickness given by observational 
uncertainty, and the intersection of plane and model 
creates a locus of possible rd and tIcI values in the host 
that could have led to the observations. Combining re- 
sults for multiple metrics further restricts the possible 
ranges of rd and nc\. This method provides good con- 
straints on host geometry from even a single epoch of 
data. Using multiple epochs of data for the same su- 
pernova, the possible host configurations can be con- 
strained even further In the future we will investigate 
other observables to find other independent constraints 
that may produce a unique rd-?ici estimate. The in- 
tersection of multiple metrics should also provide an 
indication of the precision and accuracy of our maps. 
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The results presented in this paper show the poten- 
tial for our approach in giving insight into the com- 
plicated dependence of polarization on SN geometry. 
It also provides a connection between the theoretical 
emergent Stokes spectra and robust observables, such 
as Pmax and polarized equivalent width. 

In the future, we will use this code to explore a 
larger range of parameter space to predict the obser- 
vational signature of more host configurations. Once a 
library of configurations and spectropolarimetric sig- 
natures has been developed, we will be able to under- 
stand degeneracies in the effects of a variety of host 
parameters, as well as develop methods to break these 
degeneracies with observations. Our simulations will 
also provide new observational tests for SN explosion 
codes by predicting spectropolarimetric signatures us- 
ing characteristic clump sizes and distributions that 
emerge from these models. 
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Table 2: Model parameters for investigation the effect 
of varying the number of clumps while keeping the 
clump size constant. 

Run: rci = O-lGv-max 



nci 


km s^^ 


TCF 


1 


3,200 


0.02 


4 


3,200 


0.14 


8 


3,200 


0.28 


16 


3,200 


0.61 


32 


3,200 


1.23 


64 


3,200 


2.24 


126 


3,200 


4.87 


256 


3,200 


9.76 


384 


3,200 


14.61 


512 


3,200 


19.42 


768 


3,200 


29.32 


1024 


3,200 


38.85 


Run: 


rci = 0.3^ 


^'^max 


"-cl 


7-cl 

kms~^ 


TCF 


1 


6,400 


0.08 


4 


6,400 


0.54 


8 


6,400 


1.11 


16 


6,400 


2.42 


32 


6,400 


4.90 


64 


6,400 


9.74 


126 


6,400 


19.46 


256 


6,400 


39.05 


384 


6,400 


58.45 


512 


6,400 


77.70 


768 


6,400 


117.30 


1024 


6,400 


155.42 



Table 3: Model parameters for investigation the effect 
of varying the clump radius while keeping the number 
of clumps constant. 

Run: Uci = 16 



«cl 


kms^^ 


TCF 


16 


200 


0.002 


16 


400 


0.009 


16 


800 


0.038 


16 


1,600 


0.151 


16 


3,200 


0.607 


16 


6,400 


2.429 


16 


12,800 


9.716 




Run: Tici - 


:32 


n-cl 


r-ci 
km s~^ 


TCF 


32 


200 


0.005 


32 


400 


0.019 


32 


800 


0.077 


32 


1,600 


0.306 


32 


3,200 


1.225 


32 


6,400 


4.900 


32 


12,800 


19.602 
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Table 4: Model parameters for investigation the effect 
of varying the number and radius of clumps to keep 
TCF constant. 

Run: TCF -0.15 



rici 


km s^^ 


TCF 


1 


6,400 


0.084 


4 


3,200 


0.136 


16 


1,600 


0.152 


64 


800 


0.152 


256 


400 


0.153 


1024 


200 


0.153 


Run 


: TCF- 


0.60 


"cl 


km s^^ 


TCF 


1 


12,800 


0.335 


4 


6,400 


0.543 


16 


3,200 


0.607 


64 


1,600 


0.609 


256 


800 


0.610 


1024 


400 


0.607 


Run: TCF ~ 


2.4 


"-cl 


kms^^ 


TCF 


1 


24,600 


2.11 


4 


12,800 


2.35 


16 


6,400 


2.42 


64 


3,200 


2.43 


256 


1,600 


2.44 


1024 


800 


2.43 
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